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Abstract. We bring together those systems of hydrodynamical type that can be written 
as geodesic equations on diffeomorphism groups or on extensions of diffeomorphism groups 
with right invariant L 2 or H 1 metrics. We present their formal derivation starting from 
Euler's equation, the first order equation satisfied by the right logarithmic derivative of 
a geodesic in Lie groups with right invariant metrics. 
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A fluid moves to get out of its own way as efficiently as possible. 

Joe Monaghan 

1 Introduction 

Some conservative systems of hydrodynamical type can be written as geodesic equations on the 
group of diffeomorphisms or the group of volume preserving diffeomorphisms of a Riemannian 
manifold, as well as on extensions of these groups. Considering right invariant L 2 or H 1 metrics 
on these infinite dimensional Lie groups, the following geodesic equations can be obtained: 
the Euler equation of motion of a perfect fluid [21 [10] , the averaged Euler equation [31j [50] , 
the equations of ideal magneto- hydrodynamics [MJ [32] , the Burgers inviscid equation [7] , the 
template matching equation [13 [55], the Korteweg-de Vries equation [33], the Camassa-Holm 
shallow water equation [8l 1381 [29] . the higher dimensional Camassa-Holm equation (also called 
EPDiff or averaged template matching equation) [20], the superconductivity equation [39], the 
equations of motion of a charged ideal fluid [57], of an ideal fluid in Yang-Mills field [13] and 
of a stratified fluid in Boussinesq approximation [HJ EE] • 

For a Lie group G with right invariant metric, the geodesic equation written for the right 
logarithmic derivative u of the geodesic is a first order equation on the Lie algebra g, called 
the Euler equation. Denoting by ad(u) T the adjoint of ad(u) with respect to the scalar product 
on g given by the metric, Euler's equation can be written as jj^u = — &d(u) T u. In this survey 
type article we do the formal derivation of all the equations of hydrodynamical type mentioned 
above, starting from this equation. 

By writing such partial differential equations as geodesic equations on diffeomorphism groups, 
there are various properties one can obtain using the Riemannian geometry of right invariant 
metrics on these diffeomorphism groups. We will not focus on them in this paper, but we list 
some of them below, with some of the references. 

For some of these equations smoothness of the geodesic spray on the group implies local well- 
posedness of the Cauchy problem as well as smooth dependence on the initial data. This applies 

*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
|http: / /www.emis.de/journals / SIGMA /symmetry2007.html | 
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for the following right invariant Riemannian metrics: L 2 metric on the group of volume preserv- 
ing diffeomorphisms [1U| . H 1 metric on the group of volume preserving diffeomorphisms on a 
boundary free manifold [50j . on a manifold with Dirichlet boundary conditions [31|, I51| and with 
Neumann or mixt boundary conditions [51U13j . H 1 metric on the group of diffeomorphisms of the 
circle [50|, [29] and on the Bott-Virasoro group [9] , and H 1 metric on the group of diffeomorphisms 
on a higher dimensional manifold [15J . 

There are also results on the sectional curvature (with information on the Lagrangian sta- 
bility) [21ll2lll8l[3g[Mlll6l[Ml[ia[MlE2lEZ], on the existence of conjugate points (37J HQ] 
and minimal geodesies [6] , on the f initeness of the diameter [52], [53] [11] , on the vanishing of 
geodesic distance |33j . as well as on the Riemannian geometry of subgroups of diffeomorphisms 
as a submanifold of the full dif feomorphism group |36[ d] [26] [5"5] . 

2 Euler's equation 

Given a regular Frechet-Lie group in the sense of Kriegb-Michor [28], and a (positive definite) 
scalar product ( , ) : q x q — > M on the Lie algebra q, we can define a right invariant metric 
on G by g x {^ifi) = ((.x -1 , nx~ 1 ) for £,77 E T X G. The energy functional of a smooth curve 
c : / = [a, b] — > G is defined by 



where 5 r denotes the right logarithmic derivative (angular velocity) on the Lie group G, i.e. 
5 r c(t) = c!{t)c{t)~ l € Q. We assume the adjoint of ad(X) with respect to ( , ) exists for all 
X E q and we denote it by ad(X) T , i.e. 



The corresponding notation in [3] is B(X,Y) = ad(Y) T X for the bilinear map B : q x q — > q. 

Theorem 1. The curve c : [a,b] — > G is a geodesic for the right invariant metric g on G if and 
only if its right logarithmic derivative u = 5 r c : [a, b] — ► g satisfies the Euler equation: 



E 




(ad(X)'y,Z) = {Y,[X,Z}} 



VX,Y,Z£ 9 . 



—u = — ad ti) u. 



dt 



(2.1) 
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Proof. We denote the given curve by c$ and its logarithmic derivative by uq. For any variation 
with fixed endpoints c(t,s) G G, t G [a,b], s G (—s,e) of the given curve Co, we define u = 
(c\c)c _1 and v = (d 3 c)c . In particular u(-,0) = uq, and we denote v(-,0) by vq. 
Following [35] we show first that 



dtv — d s u = [u,v]. (2.2) 

For each h G G we consider the map Fh(t,s) = (t, s,c(t, s)h) for t G [a, b] and s G (— e,e). The 
bracket of the following two vector fields on [a, b] x (— e, e) x G vanishes: 

(t, s, g) i-> d t + u(t, s)g, (t, s, g) i-> d s + v(t, s)g. 

The reason is they correspond under the mappings F^, h G G, to the vector fields 5f and <9 S on 
[a, 6] x (— e, e) (with vanishing bracket). Hence = [dt + u^r, 9 S + u^] = (dtv)g — (d s u)g — [u, v]g, 
because the bracket of right invariant vector fields corresponds to the opposite bracket on the 
Lie algebra q, so the claim (I2.2p follows. 

As in [M] we compute the derivative of E(c) = % J ( u > u)dt with respect to s, using the fact 
that v(a, s) = v(b, s) = 0. 

d s E(c) = / (d s u,u)dt = / (dtv — [u, v], u)dt = — / {v, dtu + &d(u) T u)dt. 

J a J a J a 



The curve cq in G is a geodesic if and only if this derivative vanishes at s = for all variations c 

777 ' 



of cq, hence for all : [a, b] —> g. This is equivalent to Auq = — ad(uo) T ^o- 



The Euler equation for a left invariant metric on a Lie group is = ad(u) T u. In the case 
G = 50(3) one obtains the equations of the rigid body. 

Denoting by ( , ) the pairing between g* and g, the inertia operator [3] is defined by 

A: 0^0*, A(X) = (X,-), i.e. (A(X), Y) = (X, Y), VX,Yeg. 

It is injective, but not necessarily surjective for infinite dimensional 0. The image of A is called 
the regular part of the dual and is denoted by 0* cg . 

Let ad* be the coadjoint action of on 0* given by (ad*(X)m,Y) = (m, — ad(X)Y), for 
m G 0*. The inertia operator relates ad(X) T to the opposite of the coadjoint action of X, i.e. 

ad*(X)A(Y) = -A(ad(X) T Y). (2.3) 

Hence the inertia operator transforms the Euler equation (|2.ip into an equation for m = A(u): 

— m = ad*(u)m, (2-4) 

result known also as the second Euler theorem. 

First Euler theorem states that the solution of ([2.4]) with m(a) = itiq is 

m(t) = Ad*(c(t))m , 

where u = 5 r c and c(a) = e. Indeed, = &d*(5 r c) Ad*(c)mo = &d*(u)m. 

Remark 1. Equation (|2.4|) is a Hamiltonian equation on 0* with the canonical Poisson bracket 

{fMm) = (,„,[ V*]), /,, 60 -(,-) 

and the Hamiltonian function h G C°°(0*), h(m) = ^(mjA" 1 ??!) = i(m, u). 
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Remark 2. The Euler-Lagrange equation for a right invariant Lagrangian L : TG — > M with 
value I : g — > M at the identity is: 

d 51 _ ^*t u \ M 
dt 5u 5u ' 

also called the right Euler-Poincare equation [471 130j . The Hamiltonian form ([2.4)1 of Euler's 
equation is obtained for l(u) = \{u,u) since the functional derivative is A(u) in this case. 



3 Ideal hydrodynamics 

Let G = Diff M (M) be the regular Frechet Lie group of volume preserving diffeomorphisms 
of a compact Riemannian manifold (M, g) with induced volume form ^. Its Lie algebra is 
g = X M (M), the Lie algebra of divergence free vector fields, with Lie bracket the opposite of the 
usual bracket of vector fields ad(X)Y = — [X, Y]. We consider the right invariant metric on G 
given by the L 2 scalar product on vector fields 

(X,Y) = [ g(X,Y)n. (3.1) 
JM 

In the L 2 orthogonal decomposition 3t(M) = 3£ M (M) © grad(C°°(M)), we denote by P the 
projection on £ M (M). The adjoint of ad(X) is ad(X) T y = P(V X Y + (VI) T 7) where V 
denotes the Levi-Civita covariant derivative. Indeed, 

(ad(X) T y, Z)= f g(Y, [Z, X])fi = [ g(Y, V Z X - V X Z)^ 

JM JM 
= [ g({X7X) T Y,Z)n+ [ g(V x Y,Z)n=(P(V x Y + (VX) T Y),Z), 

JM JM 

with (VX) T denoting the adjoint of the (l,l)-tensor VX relative to the metric: g(V zX,Y) = 
g(Z,(VX) T Y). In particular ad(X) T X = P(V X X) = V X X + gradp, with p the smooth 
function uniquely defined up to a constant by Ap = div(Vx^O- Now Theorem [1] assures that 
the geodesic equation in Diff^(M), in terms of the right logarithmic derivative u of the geodesic, 
is Euler's equation for ideal flow with velocity u and pressure p |4H [2l [TO]: 

dtu = — V u u — gradp, divn = 0. (3-2) 



The geodesic equation (|3.2j) written for the vorticity 2- form u = du , b denoting the inverse 
of the Riemannian lift [J and L the Lie derivative, is 

d t oj = -L u uj, (3.3) 
because (V u u)^ = L u ^ — hd(g(u,u)) and (gradp) b = dp. 



4 Burgers equation 

Let G = Diff(S' 1 ) be the group of orientation preserving diffeomorphisms of the circle and 
q = X(S' 1 ) the Lie algebra of vector fields. The Lie bracket is [X, Y] = X'Y — XY' , the negative 
of the usual bracket on vector fields (vector fields on the circle are identified here with their 
coefficient functions in C°°(S 1 )). We consider the right invariant metric on G given by the L 2 
scalar product (X,Y) = f gl XYdx on q. The adjoint of ad(Jf) is ad(X) T Y = 2X'Y + XY', 
because: 

(ad(X) T Y,Z) = I Y{X'Z -XZ')dx= f (X'Y + (XY)')Zdx = {2X'Y + XY', Z). 
Js 1 Js 1 
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It follows from Theorem[T]that the geodesic equation on Diff (S 1 ) in terms of the right logarithmic 
derivative u : I — > C co (S 1 ) is Burgers inviscid equation [7j: 

d t u = -3uu'. (4.1) 

The higher dimensional Burgers equation is the template matching equation, used for 
comparing images via a deformation induced distance. It is the geodesic equation on Diff(M), 
the diffeomorphism group of a compact Riemannian manifold (M,g), for the right invariant 1? 
metric [151155] : 

dtu = —V u u — (div u)u — i grad g(u, u). (4-2) 
Indeed, 

ad(X) T = (divX)l + V x + (VX) T , Vl€ 3£(M), (4.3) 

because as m Section [3] we compute (ad(X) T Y,Z) = f M g{(VX) T Y,Z)n + f M g(V x Y,Z)fi- 
f M L x g(Y, Z)n = {{VX) T Y + V X Y + (divX)F, Z) for all vector fields X, Y, Z on M. 

In particular for M = S 1 and u a curve in X(S' 1 ), identified with C co {S l ), divu = v! and 
g(u,u) = u 2 , so each of the three terms in the right hand side of (|4.2p is —uu' and we recover 
Burgers equation (14. ip . 



5 Abelian extensions 

A bilinear skew-symmetric map lo : q x q — > V is a 2-cocycle on the Lie algebra g with values in 
the g-module V if it satisfies the condition 

^ / u([X 1 ,X 2 ],X 3 ) = J2b(X 1 )uj(X 2 ,X 3 ), X 1 ,X 2 ,X 3 e Q , 

cycl cycl 

where b : q — > L(V") denotes the Lie algebra action on V. It determines an Abelian Lie algebra 
extension g :=V xi^g ofg by the g-module V with Lie bracket 

[(ui,Xi), (^ 2; X 2 )] = (6(Xi)«2 - 6(X 2 )^ + wpTi, X 2 ), [Xi, A 2 ]). (5.1) 

There is a 1-1 correspondence between the second Lie algebra cohomology group H 2 (q, V) and 
equivalence classes of Abelian Lie algebra extensions 0— > V ^> Q — > fl — > 0. 

When G is infinite dimensional, the two obstructions for the integrability of such an Abelian 
Lie algebra extension to a Lie group extension of the connected Lie group G involve tt\{G) 
and tt 2 (G) [43] . The Lie algebra 2-cocycle lo is integrable if 

• the period group ILj C V (the group of spherical periods of the equivariant V-valued 
2-form on G defined by ui) is discrete and 

• the flux homomorphism F w : 7Ti(G) — > H 1 (q,V) vanishes. 

Then for any discrete subg roup r of the subspace of g-invariant elements of V with T ^ 11^ , 
there is an Abelian Lie group extension 1— >T— > G ^ G — >lofGbyT = V/T. 
There are two special cases: 

1. Semidirect product: q = V X q, obtained when lo = 0. 

An example is the semidirect product g*xG for the coadjoint G-action on g*, called the 
magnetic extension in [3J. It has the Lie algebra g* X q, a semidirect product for the 
coadjoint g-action b = ad* on g*. 
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2. Central extension: g = V x w g, obtained when b = 0. 

An example is the Virasoro algebra Kx u X(S' 1 ), a central extension of the Lie algebra of 
vector fields on the circle given by the Virasoro cocycle u(X, Y) = j sl (X'Y" — X"Y')dx. 
It has a corresponding Lie group extension of the group Diff (S 1 ) of orientation preserving 
diffeomorphisms of the circle, defined by the Bott group cocycle: 

c(ip,tfj)= f logO' o<0)dlogt//, (p,i> G Diff(5 1 ). (5.2) 
Js 1 

An example of a general Abelian Lie algebra extension is C°°(M) X(M), the Abelian 
extension of the Lie algebra of vector fields on the manifold M with the opposite bracket by 
the natural module of smooth functions on M, the Lie algebra action being b{X)f = —Lxf- 
The cocycle u : X(M) x %{M) — ► C°°(M) is given by a closed differential 2-form 77 on M. If r] 
is an integral form, then there is a principal circle bundle P over M with curvature rj. In this 
case the group of equivariant automorphisms of P is a Lie group extension integrating the Lie 
algebra cocycle u: 

1 -» C°°(M,T) -» Diff(P) T -» Diff(M) [P] -» 1. (5.3) 

Here C°°(M, T) is the gauge group of P and Diff(M)rpi is the group of diffeomorphisms of M 
preserving the bundle class [P] under pullbacks (group having the same identity component as 
Diff(M)). 



6 Geodesic equations on Abelian extensions 

Following [57] we write down the geodesic equations on an Abelian Lie group extension G of G 
with respect to the right invariant metric defined with the scalar product 

((v 1 ,X 1 ),{v 2 ,X 2 )y a = {v u V2)v + {Xi,X 2 ) B (6.1) 

on its Lie algebra g = V x w g. Here ( , ) s and ( , )y are scalar products on g and V. We 
have to assume the existence of the following maps: the adjoint ad(X) T : g — > g and the adjoint 
b(X) T : V — > V for any X G g, the linear map h : V — > L s kew(fl) taking values in the space of 
skew-adjoint operators on g, defined by 

(h(v)X 1 ,X 2 ) s = (u(X u X 2 ),v) v , 
and the bilinear map l:Vx7->g, defined by 

(l(v 1 ,v 2 ),X) g = {b{X)vi,v 2 ) v . 
The diamond operation o : V x V* — ► g in [1^ corresponds to our map Z via ( , )y. 

Proposition 1. The geodesic equation on the Abelian extension G for the right invariant metric 
defined by the scalar product ()6.ip on g, written for the right logarithmic derivative (/, u), i.e. 
for curves u in g and f inV, is 

a d(u) T u-h(f)u + l(f,f), 
b(u) T f. 



d 
dt 

±s 

dt J 
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Proof. We compute the adjoint of ad(i>,X) in V » u g w.r.t. the scalar product (|6,1D 

( a d(v 1 ,X 1 ) T (v 2 ,X 2 ), (v3,X 3 )) d = ((v 2 ,X 2 ), (b(Xi)v 3 - b(X 3 ) Vl + u(Xi,X 3 ), [X 1 ,X 3 ]) i 
= {v 2 MXi)v 3 ) v + (X 2 , [X u X 3 ]) B + (v2MXi,X 3 )) v - (v 2 ,b(X 3 ) Vl ) v 
= {(b(X 1 ) T v 2 , ad(X 1 ) T X 2 + h(v 2 )X 1 - l{ Vl ,v 2 )), (v 3 , X 3 )) t 

The result follows now from Euler's equation (|2.ip . 



Remark 3. When the scalar product on V is g-invariant, i.e. (b(X)v i, v 2 )v + { v i> b(X)v 2 )v = 0, 
then I is skew-symmetric and the geodesic equation becomes 

— u = — &d(u) T u — h(f)u, 

|/ = K-)/. 

7 Geodesic equations on semidirect products 

A special case of Proposition [H obtained for uj = 0, is: 

Corollary 1. T/ie geodesic equation on the semidirect product Lie group V xi G for the right 
invariant metric defined by the scalar product (|6.ip . written for the curve (f,u) inVxig, is 

= -&d(u) T u + l(f,f), 

1/ = -b(u) T f. 
dt J y ' 

It reduces to 

—u = — &d(u) T u, 

(XL 

when the scalar product on V is g-invariant. 
Passive scalar motion 

The geodesic equation on the semidirect product C°°(M) x Diff^(M) with L 2 right invariant 
metric, written for the right logarithmic derivative (/, C°°(M) x X M (M) models passive 

scalar motion [T7] : 

d t u = -V u u - gradp, 

d t f = -df{u). (7.1) 

In this case the L 2 scalar product on C°°(M) is 3t„(M)-invariant and we apply Corollary Q] to 
get this geodesic equation. 



8 Magnetohydrodynamics 

Let A : q — > q* be the inertia operator defined by a fixed scalar product ( , } on q. The scalar 
product on the regular dual g* cg = A(g) induced via A by this scalar product in q is again 
denoted by ( , ). Next we consider the subgroup g* cg xi G of the magnetic extension g* X G, 
with right invariant metric of type (|6.ip [56] • 
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Proposition 2. If the adjoint of &d(X) exists for any X E g, then the geodesic equation on 
the magnetic extension g* cg x G with right invariant metric, written for the curve (A(v),u) in 

Sreg *Q is 

—u = — &d(u) T u + &d(v) T v, 

—v = &d(u)v. 

Proof. We have to compute the map I : g* cg x g* cg — > g and the adjoint b(X) T : g* eg — > g* eg 
for b = ad*. We use the fact (|2.3p that the coadjoint action on the image of A comes from the 
opposite of ad(-) T . Then 1{A(Y{), A{Y 2 )) = &d{Y 2 ) T Yi because 

(l(A(Y 1 ),A(Y 2 )),X) = ( & d*{X)A{Y 1 ),A(Y 2 )) = -(^(X)^, Y 2 ) = ( a d(Y 2 ) T Yx,X). 

Also the adjoint of b(X) = <id*(X) exists and b(X) T A(Y) = -A(&d(X)Y). The result follows 
now from Corollary [U ■ 

For G = 50(3) and left invariant metric on its magnetic extension g* x G one obtains 
Kirchhoff equations for a rigid body moving in a fluid. 

Let G = Diff^(M) be the group of volume preserving diffeomorphisms on a compact mani- 
fold M and g = X^(M). The regular part g* eg of g* is naturally isomorphic to the quotient 
space Q 1 (M) / d£l° (M) of differential 1-forms modulo exact 1-forms, the pairing being ([a], X) = 
J M a(X)fj,, for a G 1 (M). More precisely A(X) is the coset [X^\ obtained via the Riemannian 
metric. 

Considering the right invariant L 2 metric on the magnetic extension g* eg X G determined 
by the 1? scalar product (|3.ip on vector fields, the geodesic equations for the time dependent 
divergence free vector fields u and B are (by Proposition [2]) 

d t u = -V u u + V bB - gradp, 
8 t B = —L U B. 

We specialize to a three dimensional manifold M. The curl of a vector field X is the vector 
field defined by the relation i C urix^ = dX^ and the cross product of two vector fields X and Y 
is the vector field defined by the relation (X x Y)^ = iyixH- A short computation gives 
(curlX x X) b = ixdX^ = L x X b - dg(X,X) = (V x ^) b - \dg{X,X), hence V X X = cmlX x 
X + g grad g(X, X). The geodesic equations above are in this case the equations of ideal 
magnetohydrodynamics with velocity u, magnetic field B and pressure p \54\ I32j: 

d t u = —V u u + curli? x B — gradp, 
8 t B = —L U B. 

Magnetic hydrodynamics with asymmetric stress tensor 

Let M be a 3-dimensional compact parallelizable Riemannian manifold with induced volume 
form ji and let G = Diff^(Af) with g = 3C^(M). Each vector field X on M can be identified 
with a smooth function in C°°(M,]R 3 ), and j(X) £ C°°(M, gZ(3,M)) denotes its Jacobian. Then 
u(X,Y) = [tr(j(X)dj(Y))] G SI 1 (M) / d£t° {M) is a Lie algebra 2-cocycle on g with values in the 
regular dual g* eg . 

Considering the L 2 scalar product on the Abelian extension g* eg x^ g, we get the following 
Euler equation [5] for time dependent divergence free vector fields u and B: 

dtu = —V u u + curlS x B + ti(j(B) grad j(u)) — gradp, 
d t B = —L U B, 

modeling magnetic hydrodynamics with asymmetric stress tensor T = j(B) o j(u). 
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9 Geodesic equations on central extensions 

When V = R is the trivial g-module, then the Lie algebra action b vanishes and we get a central 
extension Rx u j defined by the cocycle uj : g x g — > R. A consequence of Proposition [T] is: 

Corollary 2. The geodesic equation on a 1- dimensional central Lie group extension G of G 
with right invariant metric determined by the scalar product ((a,X), (b,Y))g = (X,Y) g + ab on 
its Lie algebra g = R x w g is 

— u = — &d(u) T u — ak(u), a G R, 
dt 

where u is a curve in g and k E L s ^ cw (g) *s defined by the Lie algebra cocycle ui via 

{k(X),Y) =u(X,Y), VX,Yeg. 

Proof. The central extension is a particular case of an Abelian extension, so Proposition [1] 
can be applied. The linear map h : R — > L s \ iCW (g) has the form h(a)X = ak(X), because 
(h(a)Xi,X 2 ) s = au)(Xi,X2) = {ak(Xi), X 2 ) & . The g-module R being trivial, j^a = 0, so a G R 
is constant. ■ 



KdV equation 

The geodesic equation on the Bott-Virasoro group (|5.2p for the right invariant L 2 metric is the 
Korteweg de Vries equation [H]. In this case the Lie algebra is the central extension of g = 
XiS 1 ) (identified with C°°(S' 1 )) given by the Virasoro cocycle u(X,Y) = f s i(X'Y" - X"Y')dx. 
The computation w(X,Y) = -2 f gl X"Y'dx = 2 J sl X"'Ydx = (X"',Y) implies k(X) = 2X'" 
and by Corollary [2] the geodesic equation for u : I — ► C°°(S 1 ) is the KdV equation: 

d t u = -3uu' - 2au"', a G R. 



10 Superconductivity equation 

Given a compact manifold M with volume form (i, each closed 2-form r/ on M defines a Lich- 
nerowicz 2-cocycle u> v on the Lie algebra of divergence free vector fields, 

u v (X,Y)= [ V (X,Y)v. 
Jm 

The kernel of the flux homomorphism 

flux^ : X G 3£ M (M) i-» [i X fA G iJ n_1 (M, R) 

is the Lie algebra X^ X (M) of exact divergence free vector fields. On a 2-dimensional manifold 
it consists of vector fields X possessing stream functions / G C°°(M), i.e. ixfJ- = df (X is the 
Hamiltonian vector field with Hamiltonian function /). On a 3-dimensional manifold it consists 
of vector fields X possessing vector potentials A G 3£(M), i.e. ix\i = dA b (X is the curl of A). 

The Lie algebra homomorphism flux^ integrates to the flux homomorphism (due to Thurston) 
Flux /t on the identity component of the group of volume preserving diffeomorphisms: 

Flux M :Diff / ,(M) ^ J ff n - 1 (M,R)/r, Flux^) = f \i 5 r (p u ) pi]dt mod T, 

Jo 
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where ip(t) is any volume preserving diffeotopy from the identity on M to ip and T a discrete 
subgroup of if n-1 (M, M). The kernel of Flux^ is, by definition, the Lie group Diff^ x (M) of exact 
volume preserving diffeomorphisms. It coincides with Diff^(M)o if and only if H n ~ 1 (M, M) = 0. 

For 7] integral, the Lichnerowicz cocycle is integrable to Diff° x (M) [23]. When M is 3- 
dimensional, there exists a vector field B on M defined with rj = —ig/j,. The 2-form rj is closed 
if and only if B is divergence free. The integrality condition of r\ expresses as J S (B ■ n)da G Z 
on every closed surface S C M. 

The superconductivity equation models the motion of a high density electronic gas in 
a magnetic field B with velocity u: 

d t u = — V u u — au x B - gradp, a£l, (10.1) 

It is the geodesic equation on a central extension of the group of volume preserving diffeomorphisms 
for the right invariant I? metric |6H 157] . when M is simply connected. 
Indeed, 

u v (X,Y)= [ v (X,Y)n = - [ n(B,X,Y)fi= f g(X x B,Y)fi = (P(X x B),Y) 
Jm Jm Jm 

hence the map k E £ s kcw(0) determined by the Lichnerowicz cocycle lo v is k(X) = P(X x B), 
with P denoting the orthogonal projection on the space of divergence free vector fields. Now 
we apply Corollary [2j 

11 Charged ideal fluid 

Let M be an n-dimensional Riemannian manifold with Levi-Civita connection V and volume 
form (i, and r] a closed integral differential two-form. Let B be an (n — 2) vector field on M 
(i.e. B e C°°(A n ~ 2 TM)) such that rj = (—l) n ~ 2 iBH is a closed two-form. The cross product of 
a vector field X with B is the vector field XxB = (ixABl^r = (ixv)K tt denoting the Riemannian 
lift. When M is 3-dimensional, then B is a divergence free vector field with r\ = —IbH and x is 
the cross product of vector fields. 

From the integrality of rj follows the existence of a principal T-bundle tt : P — > M with 
a principal connection 1-form a on P having curvature r\. The associated Kaluza-Klein metric k 
on P, defined at a point x £ P by 

k x (X,Y) = g n{x) (T x n.X,T x 7T.Y) + a x (X)a x (Y), X,Y G T X P 

determines the volume form p, = ir*fi A a on P. 

The group Diff^(P) T of volume preserving automorphisms of the principal bundle P is an 
Abelian Lie group extension of Diff ^ (M) rp] , the group of volume preserving diffeomorphisms 
preserving the bundle class [P], by the gauge group C°°(M, T) (an extension contained in (|5.3p ). 
The corresponding Abelian Lie algebra extension 

-> C°°(M) -> X A (P) T -> X^M) -» 

is described again by the Lie algebra cocycle to : X M (M) x X M (M) — > C°°(M) given by rj. 

The Kaluza-Klein metric on P determines a right invariant I? metric on the group of volu- 
me preserving automorphisms of the principal T-bundle P. The geodesic equation written in 
terms of the right logarithmic derivative (p,u), with p a time dependent function and u a time 
dependent divergence free vector field on M, is: 

dtu = —S/ U u — pu x B — gradp, 
d t p = -dp(u). 
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It models the motion of a charged ideal fluid with velocity u, pressure p and charge density p 
in a fixed magnetic field B |57j . 

Indeed, the connection a defines a horizontal lift and identifying the pair (/, X), f £ C°°(M), 
X 6 X«(M) with the sum of the horizontal lift of X and the vertical vector field given by /, we 
get an isomorphism between the Abelian Lie algebra extension C°°(M) x^ jE^(M) and the Lie 
algebra X^P)^ of invariant divergence free vector fields on P. Under this isomorphism the I? 
metric defined by the Kaluza-Klein metric k is {(/i, Xi), (f 2 , X2)} = f M (g(Xi, X2) + fifi)^- 
The L 2 scalar product on functions is 3L^(M) invariant, i.e. b(X) is skew-adjoint. The mapping 
h : C°°(Af) -» L skcw (je M (M)) is h(f)X = P{fX x B) because: 

(h(f)X,Y) = (rj(X,Y),f)= I f(i xV )(Y)p= [ fg(XxB,Y)n=(P(fXxB),Y), 

where P denotes the orthogonal projection on the space of divergence free vector fields on M. 
The result follows from Remark [3l knowing that &d(X) T X = P(VxX). 

12 Geodesies on general extensions 

A general extension of Lie algebras is an exact sequence of Lie algebras 

O^h^fl^g^O. (12.1) 

A section s : q — > g (i.e. a right inverse to the projection q — ► g) induces the following map- 
pings PQ: 

b: fl^Der(h), b(X)f = [s(X),f], 

uj: x ^h, uj(X 1 ,X 2 ) = [s(X 1 ),s(X 2 )]-s([X 1 ,X 2 }) 

with properties: 

[b(X 1 ),b(X 2 )]-b([X 1 ,X 2 ]) = a d(u J (X 1 ,X 2 )), 
J2lj([X 1 ,X 2 ],X 3 ) = Y J HX 1 MX 2 ,X 3 ). 

cycl cycl 

The Lie algebra structure on the extension q, identified as a vector space with t) © q via the 
section s, can be expressed in terms of b and u: 

[(fuXtlifaXz)] = ([fiJ 2 ] + KX 1 )f 2 -b(X 2 )f 1 +u;(X 1 ,X 2 ),[X 1 ,X 2 ]). 

In particular for t) an Abelian Lie algebra this is the Lie bracket (|5.ip on an Abelian Lie algebra 
extension. 

We consider scalar products ( , ) g on g and ( , )(, on rj and, as in Section [HJ we impose the 
existence of several maps: ad(X) T : g — > g for any X € g, ad(/) T : rj — > f) for any / G f), 
6(A") T : fj — > fj for any X € g, as well as the linear map h : f) — > i s kcw(g) defined by 

(h(f)Xi,X 2 ) g = (uj(Xi,X 2 ), /)f,, 

and the bilinear map / : f) x () — > g, defined by 

{l(fi,f2),X) g = {b(X)f 1 J 2 ) t) . 

A result similar to Proposition Q] is: 
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Proposition 3. The geodesic equation on the Lie group extension G of G by H integra- 
ting (|12.1|) . with right invariant metric determined by the scalar product 

((h,x 1 ),(f 2 ,x 2 )) i = (/ l5 / 2 ) 6 + (x 1 ,x 2 ) g , 

written in terms of the right logarithmic derivative (p, u) is: 

&d(u) T u — h{p)u + l(p, p), 
ad(/?) T p — b(u) T p. 

13 Ideal fluid in a fixed Yang— Mills field 

Let 7r : P — ► M be a principal G-bundle with principal action a : P x G — > P and let AdP = 
P X G be its adjoint bundle. The space Q k (M, AdP) of differential forms with values in AdP 
is identified with the space Q^ or (P,g) G of G-equivariant horizontal forms on P. In particular 
C°°(M, Ad P) = C°°(P,g) G . 

We consider a principal connection 1-form a 6 ^{P, q) g on P. Its curvature r\ = da+i[a, a]. 
is an equivariant horizontal 2-form r] £ Q^ or ,(P, q) g , hence it can be viewed as a 2-form on M 
with values in Ad P. The covariant exterior derivative on g-valued differential forms on P is 
d a = X* d, with x : %(P) ~> X(P) denoting the horizontal projection, and it induces a map 
d a : Q k (M,AdP) -» O fc+1 (M,AdP). 

Let g' be a Riemannian metric on M and 7 a G-invariant scalar product on q. These data, 
together with the connection a, define a Kaluza-Klein metric on P: 

k x (X,Y) = g <x) (T x ir.X,T x ir.Y) + -y{ ax (X),a x (Y)), 1,7 6 T X P 

The canonically induced volume form on P is p = tt* pAa* det 7 , where p is the canonical volume 
form on M induced by the Riemannian metric g and a* det 7 is the pullback by a : TP — > g of 
the determinant det 7 € A dim0 $j* induced by the scalar product 7 on g. 

The gauge group of the principal bundle is identified with C QC (P,G) G , the group of G- 
equivariant functions from P to G, with G acting on itself by conjugation. The group of 
automorphisms of P, i.e. the group of G-equivariant diffeomorphisms of P, is an extension of 
Diff(M)rpi, the group of diffeomorphisms of M preserving the bundle class [P], by the gauge 
group. This is the analogue of (15. 3|) for non-commutative structure group. Restricting to volume 
preserving diffeomorphisms, we get the exact sequence: 

1 C°°(P,G) G A Diff A (P) G Diff M (M) [P] ^ 1. 

On the Lie algebra level the exact sequence is 

- C°°(P,q) g ^ X U {P) G -> 3£ M (M) -> 0. 

The horizontal lift provides a linear section : X^(M) — > X^(P) G , thus identifying the pair 
(f,X) e G°°(P,0) G ©X M (M) with X = d(f) + X hor £ X-^P) G . With this identification, the I? 
metric on X^(P) G given by the Kaluza-Klein metric can be written as 

/ K(f 1 ,X 1 ),(f 2 ,X 2 ))~p= [ g(X 1 ,X 2 )p+ [ j(h,f 2 )fi. 
Jp Jm Jp 

A particular case of a result in |T4] is the fact that the geodesic equation on the group 
Diff^(P) G of volume preserving automorphisms of P with right invariant L? metric gives the 



Jt U 



Jt P 
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equations of motion of an ideal fluid moving in a fixed Yang— Mills field. Written for the 
right logarithmic derivative (p,u) : I — > C°°(P,q) g © X M (M), these are: 

d t u = -V u u - -y(p, i u rjf - grad j>, 

d tP = -d a p(u). (13.1) 

Here u denotes the Eulerian velocity, p, viewed as a time dependent section of Ad P, denotes the 
magnetic charge, rj, viewed as a 2-form on M with values in AdP, denotes the fixed Yang-Mills 
field. The scalar product 7 being G-invariant, can be viewed as a bundle metric on Ad P. 

This result follows from Proposition Indeed, in this particular case the cocycle is u = r/ 
and the Lie algebra action is b(X)f = -df.X hoT = -d a f.X, hence b(X) T f = -b(X)f, I is skew- 
symmetric and h(f)X = P(j(f,ixf])^)- Moreover &d(X) T X = PVxX with P the projection 
on divergence free vector fields and ad(/) T / = [/, /] = 0, so (|13.ip follows. 

The equations of a charged ideal fluid from Section Qj] are obtained for the structure group G 
equal to the torus T. 

14 Totally geodesic subgroups 

Let G be a Lie group with right invariant Riemannian metric. A Lie subgroup H C G is totally 
geodesic if any geodesic c : [a,b] —> G with c(a) = e and c'(a) € rj, the Lie algebra of H, stays 
in H. 

From the Euler equation (|2.ip we see that this is the case if ad(X) T A € rj for all X S f). If 
there is a geodesic in G in any direction of f), then this condition is necessary and sufficient, so we 
give the following definition: the Lie subalgebra f) is called totally geodesic in q if ad(A) T X € f) 
for all X € f>. 

Remark 4. Given two totally geodesic Lie subalgebras rj and t of the Lie algebra g, the inter- 
section f) fl I is totally geodesic in g, but also in f) and in 6. 

Ideal fluid 

The ideal fluid flow (|3.2|) on M preserves the property of having a stream function (if M two 
dimensional), resp. a vector potential (ifMthree dimensional) if and only if Diff® x (M) isatotally 
geodesic subgroup of Diff^(M) for the right invariant L 2 metric. This means P(VxX) € 3t® x (M) 
for all X G XJf (M). 

Theorem 2 (fH]). The only Riemannian manifolds M with the property that Diff° x (M) is 
a totally geodesic subgroup o/Diff /i (M) with the right invariant L 2 metric are twisted products 
M = M. k xaF of a flat torus T k = R fc /A and a connected oriented Riemannian manifold F with 
H 1 (F, R) = 0. 

In particular the ideal fluid flow on the 2-torus preserves the property of having a stream 
function [3] and the ideal fluid flow on the 3-torus preserves the property of having a vector 
potential. 

Superconductivity 

Given a compact Riemannian manifold M, from the Hodge decomposition follows that 3t„(M) = 
j£^ x (M) £h a rm(-^0- O n a flat torus the harmonic vector fields are those with all components 
constant. 

In the setting of Section [TU1 the next proposition determines when is R xi Ul| 3L C *{M) totally 
geodesic in Mx u X M (M), for M = T 3 and rj = —isP- 
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Proposition 4 ([58J). The superconductivity equation (jlO.ip on the 3-torus preserves the pro- 
perty of having a vector potential if and only if the three components of the magnetic field B are 
constant. 

Proof. Any exact divergence free vector field X on the 3-torus admits a potential 1-form a 
with ix/J- = da, hence J" T3 g(X x B, Y)fi = J T3 iyis^ A ixfJ- = Jj3 «[y,s]jU A a. Then the totally 
geodesicity condition which, in this case, says that P(X x B) is exact divergence free for all X 
exact divergence free, is equivalent to \Y, B] = for all harmonic vector fields Y. This is further 
equivalent to the fact that the three components of the magnetic field B are constant. ■ 

Passive scalar motion 

On the trivial principal T bundle P = M x T we consider the volume form \x = fiAd9. Noticing 
that i(f,x)P- = ix^dO+ffJ,, we get the Lie algebra isomorphisms £ A (MxT) T ^ C°°(M) xX^(M) 
and 3£~ X (M x T) T = Qf 3 (M) x X° X (M), where Cfi°(M) is the subspace of functions with vanishing 
integral. 

From [55] we know that the group of equivariant volume preserving diffeomorphisms is totally 
geodesic in the group of volume preserving diffeomorphisms and from Theorem [2] we know that 
the group of exact volume preserving diffeomorphisms of a torus is totally geodesic in the group of 
volume preserving diffeomorphisms, hence by Remark|4]we obtain that for M = T 2 the subgroup 
Diffjf (Af x T) T is totally geodesic in Diff A (M x T) T . This means that C§°(M) x 3c c *(M) is 
totally geodesic in C°°(M) x Xa(M) for M = T 2 . In other words equation (|7.1j) . describing 
passive scalar motion, preserves the property of having a stream function if / has zero integral 
at the initial moment. Moreover, / will have zero integral at any moment. 

15 Quasigeostrophic motion 

Given a closed 1-form a on the compact symplectic manifold (M, a), the Roger cocycle on the 
Lie algebra X^. X (M) of Hamiltonian vector fields on M is [49] 

uj a {H f ,H g ) = [ fa(H g )a n . 
JM 

Here / and g are Hamiltonian functions with zero integral for the Hamiltonian vector fields Hj 
and H g . The integr ability of the 2-cocycle oj a to a central extension of the group of Hamiltonian 
diffeomorphisms is an open problem. Partial results are given in [24j . 

For M = T 2 the cocycle uj a can be extended to a cocycle on the Lie algebra of symplectic 
vector fields 3L a (T 2 ) by io a (d x , d y ) = io a (d x , Hf) = uj a (d y , Hj) = [27]. The extendability of u a 
to X a (M) for M an arbitrary symplectic manifold is studied in [59]. To a divergence free vector 
field X on the 2-torus one can assign a smooth function tpx on the 2-torus uniquely determined 
by X through dtpx = ix& — {ixo~) and J T2 ipxo~ = 0. Here ( ) denotes the average of a 1-form 
on the torus: (adx + bdy) = ( Jj 2 aa)dx + ( Jj 2 ba)dy. In particular ipHf — / whenever / has zero 
integral. 

Proposition 5 ([60])' The Euler equation for the I? scalar product on R x Ua X CT (T 2 ) is 

d t u = -V u u- ip u aP - gradp, (15.1) 

where the function ip u is uniquely determined by u through dip u = i u a — (i u a) and J T2 ip u o~ = 0. 

Proof. To apply Corollary [2] we compute the map k corresponding to the cocycle u a . Using 
the fact that oj a {d x ,X) = u a (d y ,X) = for all X G 3c a (T 2 ), we get 

u> a (u,X) =u a (H i!u ,X) = [ i> u a{X)a= I g( 1 p u aKx)a = (P(tl; u J),X), 

JT 2 JT 2 
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hence k(u) = P(ijj u a^). Knowing also that &d(u) T u = P(V u u), we get (|15.ip as the Euler 
equation for a = 1. ■ 

Proposition 6 Q60J). // the two components of the 1-form a on T 2 are constant, then equa- 
tion (|15.ip preserves the property of having a stream function, i.e. R x UJa X° X (T 2 ) is totally 
qeodesic in R x„ , X (T (T 2 ). In this case the restriction of (|15,ip to Hamiltonian vector fields is 

d t H^ = -V ff /^ - V« S - gradp. (15.2) 

Proof. By Theorem [2] on the 2-torus P(VjjfX) is Hamiltonian for X Hamiltonian, hence the 
totally geodesicity condition in this case is equivalent to the fact that Pfyxofi) is Hamiltonian 
for X Hamiltonian. By Hodge decomposition this means tpx^ is orthogonal to the space of 
harmonic vector fields, so 

(P(i) X c$),Y)= ( g(ip x a i ,Y)o-= [ a(Y)^ x o- = 0, V Y harmonic. 
Jt 2 Jt 2 

On the torus the harmonic vector fields Y are the vector fields with constant components and 
the functions tpx have vanishing integral by definition, so the expression above vanishes for all 
constant vector fields Y if the 1-form a has constant coefficients. ■ 

On the 2-torus with o = dx Ady and u = H^, the vorticity 2-form is du^ = d{H^f = (Aip)a, 
hence u = is the vorticity function. Since L u (du ) = LhAoju) = {Lh^oj)^ = {u),ifi}o~, the 
vorticity equation (13.3|) written for the vorticity function uj becomes 

d t oj = -{uj,ip}. 

For a = f3dy, (3 E R, we have d(ipa^)^ = dif) A a = (pd x ip)a. Hence the Euler equation (|15.2p 
written for the vorticity function u = Aip with ip the stream function of u, is the equation for 
quasigeostrophic motion in /3-plane approximation [62j EI] 

d t uj = ~{u,ip} - Pd x tp, 
with (3 the gradient of the Coriolis parameter. 



16 Central extensions of semidirect products 

Let g be a Lie algebra with scalar product ( , } g and V a g-module with g-action b and g-invariant 
scalar product ( , }y. Each Lie algebra 1-cocycle a G Z l {$,V) (i.e. a linear map a : g — > V 
which satisfies a([Xi,X2]) = b{Xi)a{X2) — b{X2)a(X\)) defines a 2-cocycle uj on the semidirect 
product V x g [35]: 

co((v 1 ,X 1 ),(v 2 ,X 2 )) = (a(X 1 ),v 2 ) v - {a(X 2 ),v 1 ) v . (16.1) 

Proposition 7. The Euler equation on the central extension (g tx V) x u R with respect to the 
scalar product ( , ) + ( , )y, written for curves u in g and f in V , is 

— &d(u) T u + aa T (f), 
b{u)f — aa{u), a 6 R, 

where a 1 : V — > g is the adjoint of a : g — > V . 



dt 
dt J 



16 



C. Vizman 



Proof. The map k € £ s kcw(^ ^ q) defined by u is k(v,X) = (a(X), — a T (v)) because 

u((vi,X 1 ),(v 2 ,X 2 )) = {a(Xi),v 2 )v ~ {a T (v 1 ),X 2 ) s = ((a(Xi), -a T (v^), (v 2 , X 2 )) Vy , s . 

The result follows from Corollaries Q] and [2j ■ 

Remark 5. More generally, a 1-cocycle a on q with values in the dual g-module V* defines 
a 2-cocycle on V x Q by 

u((v 1 ,X 1 ),(v 2 ,X 2 )) = (a(X 1 ),v 2 ) - (a(X 2 ),v 1 ), 

where ( , ) denotes the pairing between V* and V. 

17 Stratified fluid 

Let M be a compact Riemannian manifold with induced volume form [i. Let a be a closed 1-form 
on M. Then a : X(M) — > C°°(M) is a Lie algebra 1-cocycle with values in the canonical X(M)- 
module C°°(M). The 1? scalar product is £^(M)-invariant, so a defines by (|16.ip a 2-cocycle co 
on the semidirect product Lie algebra C°°{M) x X^(M): 

Lo((f 1 ,X 1 ),(f 2 ,X 2 ))= f f 2 a{X{)n- [ fia(X 2 )». (17.1) 

JM JM 

Proposition 8. The Euler equation on (C°°(M) x X M (M)) x w R with I? scalar product is 
dtu = —V u u + afcft — gradp, 

d t f = -L u f - aa(u), a € R, (17.2) 
with V the Levi-Civita covariant derivative and jj the Riemannian lift. 

Proof. We apply Proposition for g = X M (M) and F = C°°(M). In this case 6(X)/ = 
-L x f and ad(X) T X = PV X *. We compute (a T (f),X) s = J M fa(X)fi = f M g{fa$,X)fi = 
(P(fa 9 ), X) s , for all X 6 X M (M), hence a T (/) = P(/a s ). ■ 

Because <i(/a") = df A a, the equation (|17.2p written for vorticity 2-form a; = fin 1, becomes 

dfto = —L u uj + df A a, 
d t f = -L u f - aa(u). 

Proposition 9 ([SB]). Given a 2-cocycle u) determined via (|17.ip by the constant 1-form a on the 

torus M = T 2 , we have (3^ X (M) x C^(M)) x„ R is totally geodesic in (£ M (M)xC°°(M))x w R, 
where Cq°(M) is the subspace of functions with vanishing integral. 

Proof. We know from Section [H that for M = T 2 , 3t° x (M) x C^(M) is totally geodesic in 
3£fj,(M) x C°°(M). But a(u) has zero integral for u exact divergence free, a being closed. We 
have to make sure that fa^ is orthogonal to the space of harmonic vector fields for all / with 
zero integral, i.e. for all functions / such that //x is exact (ffi = du). But J M g(fa^,Y)fi = 
J M a(Y)dv = — J M Lyd A v = because Lya = for all harmonic vector fields Y on the 
2-torus, a being a constant 1-form. ■ 

Hence on the 2-torus, for constant a and initial conditions uq Hamiltonian vector field and fo 
function with zero integral, u will be Hamiltonian and / will have zero integral at every time t. 
The Hamiltonian vector field is = d y ipd x — d x ipd y and the Poisson bracket Lh 4 J = {f,ip} 
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is the Jacobian of / and tp. If a = —f3dy and a = — 1 we get the equation for stream function ip 
and vorticity function u = Atjj: 

d t uJ = - (3d x f, 

d t f = -{f,iA + Pd x ^. (17.3) 

Let £ = 9 P ~p^° be a buoyancy variable measuring the deviation of a density p from a back- 
ground value po) with g the gravity acceleration. The background stratification po is assumed 
to be exponential, characterized by the constant Brunt- Vaisala frequency N = (— g rfl ^ w )2. 
The equation for a stratified fluid in Boussinesq approximation [61] is the geodesic equa- 
tion (fTO) for /3 = N constant and / = iV" 1 ^ 

d t uj = -{u,tp} - d x £, 
dtt = -{Z^} + N 2 d x ^. 

When the Brunt -Vaisala frequency N is an integer and £ has zero integral (at time zero), 
then the stratified fluid equation is a geodesic equation on a Lie group [58] . 

18 H 1 metrics 

Camassa— Holm equation 
The Camassa Holm equation |8j 

d t (u - u") = -3uu' + 2u'u" + uv!" (18.1) 

is the geodesic equation for the right invariant metric on Diff (S 1 ) given by the H 1 scalar product 
(X,Y) = J S1 (XY + X'Y')dx = f sl X(l-d%)Ydx 
Indeed, one gets from 

(&d(X) T Y, Z) = (Y, X'Z - XZ') = [ (Y(X'Z - XZ') + Y'(X"Z - XZ"))dx 

Js 1 

= I Z(2YX' + Y'X -2Y"X' -Y"'X)dx 
Js 1 

that &d(X) T Y = (1 - dl)~ l {2YX' + Y'X - 2Y"X' - Y"'X). Plugging &&(X) T X = (1 - 
d%)- 1 (3XX'-2X'X"-XX"') into Euler's equation (J57Q) one obtains the Camassa-Holm shallow 
water equation for u : I — > C co (S 1 ). 

Since m = A(u) = u — u" , the Hamiltonian form of the Camassa-Holm equation is 

dtm = —urn! — 2u'm. 

Remark 6. Considering the right invariant H 1 metric on the Bott-Virasoro group (|5.2p . an 
extended Camassa-Holm equation is obtained [3H] 

d t (u - u") = -3uu' + 2u'u" + uv!" - 2au"', a € R. 

Indeed, the identity u(X, Y) = {k(X),Y} for the Virasoro cocycle u(X, Y) = 2 J gl X"'Ydx and 
the H 1 scalar product implies k(X) = 2(l — d^)~ 1 X'". Now by Corollary [2] the geodesic equation 
is the extended Camassa-Holm equation above. 

The homogeneous manifold Diff (S 1 )/S 1 is a coadjoint orbit of the Bott-Virasoro group. The 
Hunter Saxton equation describing weakly nonlinear unidirectional waves [22] 

d t u" = -2u'u" - uu'" 

is a geodesic equation on Diff(S' 1 )/S' 1 with the right invariant metric defined by the scalar 
product (X,Y) = J sl X'Y'dx fk\. 
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Higher dimensional Camassa— Holm equation 

The higher dimensional Camassa— Holm equation (also called EPDiff or averaged template 
matching equation) |19l [20] is the geodesic equation for the right invariant H 1 metric 

(X,Y)= [ (g(X,Y) + a 2 g(VX,VY))^ (18.2) 
JM 

on Diff(M) for compact M. Because V*V = A + Ric, this scalar product can be rewritten 
with the help of the rough Laplacian Ar = A + Ric as (X,Y) = f M (g(X — a 2 A^X, Y)[i, 
so the momentum density of the fluid m = A(u) is m = (1 — o?Ar)u. It follows that the 
adjoint of &d(X) with respect to (|18.2|) is conjugate by 1 — ot 2 An to the adjoint of ad(X) 
with respect to the I? metric (|3.ip computed to be (|4.3[) . Hence (1 — a 2 Ar) ad(X) T Y = 
(VX) T (F - a 2 A R Y) + V X (Y - a 2 A R Y) + (divX)(Y - a 2 A R Y). We get as geodesic equation 
the higher dimensional Camassa-Holm equation 

dt(l — o. 2 Ar)u = —udivu + a 2 (divu)ARU — V u u + a 2 V u (Aru) 
- {Vu) T u + a 2 (Vu) T Aru. 

In Hamiltonian form this equation is 

dtm = —V u m — (S/u) T m — (divn)m for m = A{u) = u — 0( 2 Aru. 

In particular for M = S 1 we get the Camassa-Holm equation (|18.ip . 

When M is a manifold with boundary and we put Neumann or mixed conditions on the 
boundary, then the H 1 scalar product (|18.2|) has to be replaced by 

(X,Y)= [ (g(X,Y) + 2a 2 g(DeiX,DefY)) f i, (18.3) 

JM 

where Def X = \{VX + (VI) T ) denotes the deformation (1, l)-tensor of X [15]. 
Averaged Euler equation 

For a compact Riemannian manifold M we consider the right invariant metric on the group 
Diffu(M) of volume preserving diffeomorphisms given by the H 1 scalar product (|18.2p on vector 
fields. The geodesic equation is the (Lagrangian) averaged Euler equation p3T|l50]. also called 
LAE-a equation: 

d t (l - a 2 A R )u = -V u (l - a 2 A R )u + a 2 (X7u) T (A R u) - gradp. (18.4) 

Indeed, from (ad(X) T y, Z) = f M g(Y - a 2 A R Y, V Z X - V x Z)^i = J M g((VX) T (Y - a 2 A R Y) + 
V X {Y - a 2 A R Y), Z)fjL we obtain that 

(1 - a 2 A R )( a d(X) T Y) = P((VX) T Y + V X (1- a 2 A R )Y - a 2 (VX) T (ArY)) 

and we use Euler 's equation (|2.ip to get (|18.4p . In Hamiltonian form this equation is 

dtm = —V u m — (Vu) T m — gradp 

for m = A(u) = u — cx 2 Aru. 

As in the higher dimensional Camassa-Holm equation, when Neumann or mixed conditions 
on the boundary of M are imposed, one has to consider the H 1 scalar product (|18.3p . 
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19 Systems of two evolutionary equations 

From [T2] we know that a basis for # 2 (£(S' 1 ), C 00 ^ 1 )) is represented by a(X,Y) = X'Y - 
XY' and the Virasoro cocycle u(X,Y) = f s i(X'Y" - X"Y')dx € R C C°°{S l ); a basis for 
H 2 (X(S 1 ),H 1 (S 1 )) is represented by the cocycles 

o-i (X, Y) = XY" - X"Y, cji(X, Y) = X'Y" - X"Y'; 

a basis for H 2 (X(S 1 ), 2 (5 1 )) is represented by the cocycles 

a 2 {X, Y) = X"'Y - XY'", u 2 {X, Y) = X"'Y' - X'Y'". 

Only the cocycles u, u\ and oj 2 (whose expressions involve only derivatives of X and Y) integrate 
to group cocycles [15] . 

The Euler equations for the 1? or H 1 scalar product on the corresponding Abelian extensions 
provide systems of two equations, generalizing Burgers (|4,ip or Camassa-Holm ()18.ip equation. 
We exemplify with the 2-cocycle a taking values in the module of functions on the circle. The 
Euler equations for the L 2 scalar product on C ao (S 1 ) xi X(S' 1 ) and on C°°(S 1 ) » g X(S 1 ) are 

d t u = -3uu' - ff, 
dtf = ~uf - u'f, 

and 

d t u = - W + uf + 2u'f - ff, 
dtf = ~uf - u'f. 

The Euler equations for the H 1 scalar product on C°°(S' 1 ) xi X(S' 1 ) and on C co (S 1 ) xi CT X(S' 1 ) 
are 

d t ( u - u") = -3uu' + 2u'u" + uu'" - ff + ff", 
d t (f-f") = -uf'-u'f + uf"' + u'f", 

and 

d t {u - u") = -3uu' + 2u'u" + uu'" - 2u'f - uf + 2u' f" + uf" - ff + ff", 
d t (f-f") = -uf'-u'f + uf"' + u'f". 

One can consider central extensions of semidirect products of j£(S' 1 ) with modules of densities 
as in Remark[5] j3S]- For instance the 1-cocycle a{X) = X" on X(S 1 ) with values in the 
module dual to C°°(S 1 ), gives the 2-cocycle w((/ l9 JTi), (f 2 ,X 2 )) = J s i(X'{f 2 - X%fi)dx on the 
semidirect product C co (S 1 ) x X(S' 1 ). The geodesic equation for the L 2 scalar product on the 
central extension (C 00 (S' 1 ) x X(S' 1 )) x w lis 

d t u = -3uu' - ff - af", 

d t f = -uf - u'f - au", a £ R. 

20 Conclusions 

This survey article presents the formal deduction as geodesic equations on diffeomorphism groups 
with right invariant metrics of several PDE's of hydrodynamical type. Sometimes extensions of 
diffeomorphism groups by central or Abelian sugroups come into play and the corresponding 
Lie algebra 2-cocycles introduce additional terms to the geodesic equations. 

These equations are Hamiltonian equations too, possessing rich geometric structures. Some 
of them are completely integrable. But presenting these results is beyond the scope of this 
article. 
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